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ABSTRACT 

Echelle spectrophotometry of the planetary nebula NGC 5307 is presented. 
The data consists of VLT UVES observations in the 3100 to 10360 A range. 
Electron temperatures and densities have been determined using different line 
intensity ratios. We determine the H, He, C, and O abundances based on re- 
combination lines, these abundances are almost independent of the temperature 
structure of the nebula. We also determine the N, O, Ne, S, Cl, and Ar abun- 
dances based on collisionally excited lines, the ratios of these abundances relative 
to that of H depend strongly on the temperature structure of the nebula. From 
the O II/[0 III] ratios we find a t 2 = 0.056 ± 0.005. The chemical composition of 
NGC 5307 is compared with those of the Sun and the Orion nebula. From the 
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study of the relative intensities of the O II recombination lines of multiplet 1 in 
this and other nebulae it is found that for electron densities smaller than about 
5000 cm" 3 collisional redistribution is not complete, this effect has to be taken 
into account to derive the O abundances for those cases in which not all the lines 
of the multiplet are observed. From the A 4649 O II versus iV e (Cl III) diagram 
we find a critical electron density of 1325 cm" 3 for collisional redistribution of 
the O II lines of multiplet 1. Also based on this diagram we argue that the O II 
and the [O III] lines originate in the same regions. We also find that the radial 
velocities and the FWHM of the O II and [O III] lines in NGC 5307 are simi- 
lar supporting the previous result. These two results imply that for NGC 5307 
and probably for many other gaseous nebulae chemical inhomogeneities are not 
responsible for the large temperature fluctuations observed. 

Subject headings: planetary nebulae: abundances — planetary nebulae: individual 
(NGC 5307) 

1. Introduction 

The main aim of this paper is to make a new determination of the chemical abun- 
dances of NGC 5307 including the following improvements over previous determinations: 
the consideration of the temperature structure that affects the helium and heavy elements 
abundance determinations, the derivation of the O and C abundances from recombination 
line intensities, the consideration of the collisional excitation of the triplet He I lines from 
the 2 3 S level by determining the electron density from many line intensity ratios, and the 
study of the 2 3 S level optical depth effects on the intensity of the triplet lines by observing a 
large number of singlet and triplet lines of He I. To determine the temperature structure it is 
also important to discuss if the O II and the [O III] lines originate in the same regions under 
the same physical conditions. To this effect we will study the relative line intensities of the 
O II multiplet 1 in this object and other gaseous nebulae, moreover for the same reason we 
will also consider the wavelengths and line widths of the O II and [O III] lines in NGC 5307. 

In sections 2 and 3 the observations and the reduction procedure are described. In 
section 4 temperatures and densities are derived from five and four different intensity ratios 
respectively; also in this section, the mean square temperature fluctuation, t 2 , is determined 



1 Bascd on observations collected at the European Southern Observatory, Chile, proposal number ESO 
68.C-0149(A). 
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from the O II/ [O III] line intensity ratios. In section 5 ionic abundances are determined 
based on recombination lines that are almost independent of the temperature structure, and 
ionic abundances based on ratios of collisionally excited lines to recombination lines that do 
depend on the temperature structure of the nebula. In section 6 the total abundances are 
determined and compared with other determinations for the same object. In sections 7 and 
8 we present the discussion and the conclusions. 

2. Observations 

The observations were obtained with the Ultraviolet Visual Echelle Spectrograph, UVES 
(D'Odorico et al. 2000), at the VLT Kueyen Telescope in Chile. We observed simultaneously 
with the red and blue arms, each in two settings, covering the region from 3100 A to 10360 
A (see Table 1). The wavelength regions 5783-5830 AA and 8540-8650 AA were not observed 
due to the separation between the two CCDs used in the red arm. There were also two small 
gaps that were not observed, 10084-10088 AA and 10252-10259 AA, because the two redmost 
orders did not entirely fit within the CCD. In addition to the long exposure spectra we took 
60 second exposures for the four observed wavelength ranges to check for possible saturation 
effects. 

The slit was oriented east-west and the atmospheric dispersion corrector (ADC) was 
used to keep the same observed region within the slit regardless of the air mass value. The 
slit width was set to 3.0" and the slit length was set to 10" in the blue arm and to 12" 
in the red arm; the slit width was chosen to maximize the S/N ratio of the emission lines 
and to maintain the required resolution to separate most of the weak lines needed for this 
project. The FWHM resolution for the NGC 5307 lines at a given wavelength is given by 
A A ~ A/8800. The reductions were made for an area of 3" x 10" . The center of the slit 
was placed 5" south of the central star, a region of high emission measure. The dimensions 
of the object are 18.8" x 12.9" (e. g. Tylenda et al. 2003). The average seeing during the 
observations amounted to 0.8". 

The spectra were reduced using the IRAF 2 echelle reduction package, following the 
standard procedure of bias subtraction, aperture extraction, flatfielding, wavelength calibra- 
tion and flux calibration. For flux calibration the standard star EG 247 was observed; the 
observing time amounted to 600 seconds in each of the four observed wavelength ranges. 



2 IRAF is distributed by NOAO, which is operated by AURA, under cooperative agreement with NSF. 
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3. Line Intensities, Reddening Correction and Radial Velocities 

Line intensities were measured integrating all the flux in the line between two given 
limits and over a local continuum estimated by eye. In the few cases of line-blending, the 
line flux of each individual line was derived from a multiple Voigt profile fit procedure. All 
these measurements were carried out with the splot task of the IRAF package. 

The reddening coefficient, C(H/3) = 0.59 ±0.10 dex, was determined by fitting the 
observed 7(H/3)/7(H Balmer lines) ratios to the theoretical ones computed by Storey & 
Hummer (1995) for T e = 10,000 K and N e = 2000 cm" 3 , (see below) and assuming the 
extinction law of Seaton (1979). 

Table 2 presents the emission line intensities of NGC 5307. The first and second columns 
include the observed wavelength in the heliocentric framework, A, and the adopted laboratory 
wavelength, A . The third and fourth columns include the ion and the multiplet number, or 
the Balmer, B, or Paschen, P, transitions for each line. The main sources used for the line 
identifications and laboratory wavelenghts were: Esteban et al. (1998); Hyung et al. (2000); 
Liu et al. (2001); Wiese et al. (1996); Pequignot & Baluteau (1988), and references therein. 
The fifth and sixth columns include the observed flux relative to H/3, -F(A), and the flux 
corrected for reddening relative to to H/3, /(A). To combine all the line intensities, from the 
four different instrumental settings, on the same scale, we multiplied the intensities in each 
setting by a correction factor obtained from the lines present in more than one setting; these 
corrections were all smaller than 2%. The seventh column includes the fractional error (la) 
in the line intensities; these errors were estimated from a signal-to-noise versus signal curve 
obtained from a fit to the dispersion in the observed to predicted line intensities for all the 
measured hydrogen lines. 

A total of 184 distinct emission lines were measured; of them 137 are permitted, 46 are 
forbidden and 1 is semiforbidden (see Table 2). In this tally only the main component of 
the blended lines is included. Many lines redward of 7000 A inclunding several Paschen lines 
and A 9531 [S III] were detected, but they were strongly affected by atmospheric features in 
absorption and emission rendering their intensities unreliable. 

We derived the heliocentric radial velocity of NGC 5307 based on the 14 strongest H I 
and He I lines and amounts to +29.4 ± 0.1 km sec -1 . This value is smaller than that 
derived by Meatheringham, Wood, & Faulkner (1988) that amounts to +40 ± 4 km sec -1 , 
the difference probably is due to the different areas observed by the two groups. 

The center of our observing slit was located at about half the distance between the 
central star and the outer edge of the main body of the nebula. Most of the emission lines 
show a single line profile because we are observing material moving from a tangential velocity 
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up to an angle of 60 degrees relative to the plane of the sky. The exceptions are the blue 
and red components of the AA 6313 [S III], 6716 [S II], 6731 [S II], and 9069 [S III] lines 
that show a splitting of 35 ± 3 km sec -1 , because they mainly originate in the outer regions 
of the nebula. Dividing this value by the sine of 60 degrees we find a 2V exp value of about 
40 km sec -1 of this material relative to the central star. This value is similar to the 2V exp 
values of 30 km sec -1 and 21.2 km sec -1 found by Acker (1976) and Meatheringham, Wood, 
& Faulkner (1988) respectively. 



4. Physical Conditions 

4.1. Temperatures and densities 

The temperatures and densities presented in Table 3 were derived from the line inten- 
sities presented in Table 2. Most of the determinations were carried out based on the IRAF 
subroutines. 

The contribution to the intensities of the AA 7319, 7320, 7331, and 7332 [O II] lines due 
to recombination was taken into account based on the following equation: 

7^(7319 + 7320 + 7331 + 7332)/J(H/3) = 9.36(T/10 4 ) a44 x 0++/H+, (1) 

(see Liu et al. 2000). 

The contribution to the intensity of the A 5755 [N II] line due to recombination was 
taken into account based on the following equation: 

J i? (5755)//(H/3) = 3.19(T/10 4 ) - 30 x N++/H+, (2) 

(see Liu et al. 2000). 

The Balmer continuum temperature was determined from the following equation: 

T = 368 x (1 + 0.259y + + 3.409y ++ )(Bac/Hll)- 3/2 K, (3) 
(see Liu et al. 2001) where Bac/Hll is in A -1 , y + = He+/H+, and y ++ = Hc+ + /H+. 



4.2. Temperature variations 

To derive the ionic abundance ratios the average temperature, T , and the mean square 
temperature fluctuation, t 2 , were used. These quantities are given by 

TJN N) = Jr.frWrWMriV 
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and 



e = J(T e -T ) 2 N e NidV 



T 2 / N e N t dV ' (5) 
respectively, where N e and iVj are the electron and the ion densities of the observed emission 
line and V is the observed volume (Peimbert 1967). 

To determine To and t 2 we need two different methods to derive T e : one that weighs 
preferentially the high temperature regions and one that weighs preferentially the low tem- 
perature regions. In this paper we have used the temperature derived from the ratio of the 
[O III] AA 4363, 5007 lines, T( 4 3 6 3/5 007 ), that is given by 



T 



(4363/5007) 



T 



1 + 



1 /90800 



T n 



3 \t 2 



(6) 



and the temperature derived from the ratio of the recombination lines of multiplet 1 of O II 
to the collisionally excited lines of [O III] that is given by 



T(0 n rec/O in coll) — T( 46 49/ 5 007) — Hnx). (7) 

From these equations we obtain that for 0++(R/C): t 2 = 0.056 ±0.005, and T = 10100 
K. Since about 90 % of the oxygen is twice ionized (see below) we will use this t 2 value as 
representative for the whole object. From T(Bac), T(0 II), T(0 III), and equations 6-9 by 
Peimbert et al. (2000) we obtain t 2 = 0.03 ± 0.05 this value is not accurate enough to be 
useful in the abundance determinations. 



5. Ionic Chemical Abundances 

5.1. Helium ionic abundances 

To obtain He + /H + values we need a set of effective recombination coefficients for the He 
and H lines, the contribution due to collisional excitation to the helium line intensities, and 
an estimate of the optical depth effects for the helium lines. The recombination coefficients 
used were those by Storey & Hummer (1995) for H, and those by Smits (1996) and Benjamin, 
Skillman, & Smits (1999) for He. The collisional contribution was estimated from Sawey & 
Berrington (1993) and Kingdon & Ferland (1995). The optical depth effects in the triplet 
lines were estimated from the computations by Benjamin, Skillman, & Smits (2002). 

Table 4 presents the He + /H + values obtained from the eight best observed helium 
lines for N e = 2500 ± 500. From a maximum likelihood method we obtained He + /H + = 
0.09055 ± 0.00151, r 3889 = 0.93 ± 0.57, and t 2 = 0.031 ± 0.014 
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From the 7(4686) /J(H/3) value presented in Table 2, T e = 13500 K, and the recombi- 
nation coefficients by Brocklehurst (1971) we obtain iV(He ++ )/iV(H + ) = 0.00726 ± 0.00018. 
The helium ionic abundances are presented in Table 5. 

5.2. C and O ionic abundances from recombination lines 

The C ++ abundance was derived from the A4267 A line of C II and the effective re- 
combination coefficients computed by Davey, Storey, & Kisielius (2000) for Case A and T = 
10000 K. 

The ++ abundance was derived from the lines of multiplet 1 of O II (see Figure 1) 
together with the effective recombination coefficient for the multiplet computed by Storey 
(1994) under the assumption of Case B for T e = 10000 K and N e = 2000 cm" 3 . 

Six of the eight lines of multiplet 1 were observed and are presented in Table 2, the 
AA 4674 and 4696 A lines were too weak to be observed, we computed their intensities as 
follows: a I(4674)=0.019 value was determined from 1(4651) and the ratio of their Einstein A 
coefficients, since both lines originate from the same upper level; similarly a I(4696)=0.008 
value was determined from 1(4639) + 1(4661) and the ratio of their A coefficients, since the 3 
lines originate from the same upper level. The A values were taken from Wiese et al. (1996). 
Additionally the A4641.81 A line of O II is blended with the A4641.85 A line of N III, based 
on the observed intensities and the A values of other lines of the same multiplets originating 
in the same upper energy levels we estimated that the intensities of these lines, in the units 
of Table 2, amount to 0.148 and 0.053 respectively. 

The ++ abundance is almost independent of the case assumed, the difference in the 
++ /H + value between Case A and Case B is smaller than 4%. Peimbert (2003) found that 
in 30 Doradus the O II lines of multiplet 1 are in Case B based on the observed intensities of 
multiplets 19, 2, and 28 of O II that are strongly case sensitive; Peimbert, Storey, & Torres- 
Peimbert (1993) also found that the O II lines in the Orion nebula are in Case B. The line 
intensity ratios within multiplet 1 do not follow the predictions for collisional redistribution, 
this result will be discussed in section 7. Figure 1 provides an excellent visual reference to 
estimate the quality of the data, it includes 2 lines four orders of magnitude fainter than H/3 
and also shows that lines separated by 2A are completely resolved. 

The C +3 abundance was derived from: a) A4647 A , one of the three lines of multiplet 
1 of C III, b) the effective recombination coefficient used by Liu et al. (2000) for the whole 
multiplet, and c) the assumption that 7(4647 + 4650 + 4651)/7(4647) = 1.8, the LTE ratio 
for multiplet 1 (Wiese et al. 1996). 
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Most of the intensity of each of the O III and N III lines present in Table 2 is due to the 
well known Bowen resonance fluorescence mechanism (Bowen 1934) making these lines and 
are not suitable for accurate abundance determinations. 



5.3. Ionic abundances from collisionally excited lines 



The values presented in Table 6 for t 2 = 0.00 were derived with the IRAF task abund. 
For N II, O II, and S II we used T e = 11000 K and N e = 3500 cm" 3 . For O III, Ne III, S III, 
CI III, CI IV, Ar III, and Ar IV we used T e = 11800 K and N e = 2500 cm" 3 . For Ar V and 
K IV we used T e = 13500 K and iV e = 2500 cm" 3 . 

To derive the abundances for t 2 = 0.056 we used the abundances for t 2 = 0.00 and 
the following equations for t 2 > 0.00 presented by Peimbert (1967) and Peimbert & Costero 
(1969): 
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where a is the dependence on the temperature for a given recombination line, AE CEL is the 
energy difference between the ground and the excited level, and [iV(CEL)/A^(RL)](4 36 3/ 5 oo7) 
is the abundance ratio based on the temperature obtained from J(4363)//(5007), note that 
the temperature exponents presented by Peimbert & Costero (1969) for equation (8) are 
misprinted. 

To derive abundances for other t 2 values it is possible to interpolate or to extrapolate 
from the values presented in Table 6. 



6. Total Abundances 



Table 7 presents the total gaseous abundances of NGC 5307 for t 2 = and t 2 = 0.056. 
To derive the total gaseous abundances the set of equations presented below was used, where 
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the ionization correction factors, 7CF's, correct for the unseen ionization stages. 
The total He/H value is given by: 

iV(He) _ 7V(He + ) + 7V(He ++ ) 
7V(H) ~ 7V(H+) ' 

The C abundance is given by: 



(12) 



7V(H) " ICF{C) iV(H+) • (13) 

The ICF(C) was estimated from models computed with CLOUDY (Ferland 1996; Ferland 
et al. 1998) with similar O and He ionization structure and amounts to 1.05, a value very 
close to I. Half of the correction is due to the expected C + contribution and half to the C +4 
contribution. 

The following equation has been used often to obtain the N/H abundance ratio (Peim- 
bert & Costero 1969): 

7V(H) 1 ; iV(H+) 7V(0+)7V(H+)' 1 ' 

From our observations of NGC 5307, JCF(N) is in the 30 to 40 range, a very large and 
consequently uncertain correction. Relano et al. (2002); Luridiana et al. (2002) have obtained 
ICF(N + ) values bigger than those given by equation (14). Therefore we decided to use a 
different equation to obtain N/H given by: 

7V(N) _ / N(Q) \ / iV(N++) + iV(N+ 3 ) \ 

7V(H) {n(E+)) v {n(0++) + N(0^)) uv > ^ > 

where the first term on the right hand side was obtained from Table 7 and the N and O ionic 
abundances were obtained from: a) the line intensities of AA 1401 O IV, 1486 N IV], 1660 
+ 1666 O III], and 1750 N III] presented by Kingsburgh & Barlow (1994), b) the atomic 
parameters presented in the compilation by Mendoza (1983) for AA 1660 + 1666 O III], and 
the atomic parameters in IRAF for the other lines, and c) a T e of 12500 K. The values derived 
from equation (15) are those presented in Table 7. 

We prefer equation (15) to equation (14) because about 97 - 99% of the N an O atoms 
are expected to be doubly or triply ionized, while only about 1 - 3 % of the N and O atoms 
are expected to be singly ionized. 



It is also possible to determine N/H from: 

N(N) _ [7V(N++) + N{N +3 )} UV 
7V(H) " [N(H+)] V ■ 



(16) 



we prefer equation (15) to equation (16) because there is a mismatch between the abundances 
derived from UV /V line intensity ratios and those derived from the UV/UV or V/V ratios. 
For example, from the UV data by Kingsburgh & Barlow (1994) and their I(Hj3) value 
we obtain N(0) uv / N(H) V = 12.65 x 10~ 4 for t 2 = 0.00, while from their visual data we 
obtain N(0) v /N(R) v = 4.3 x 10~ 4 , also for t 2 = 0.00, a difference of a factor of three. This 
difference could be due to errors in the C(H/5) value, errors in the adopted temperature, and 
differences in the slit size and position between the UV and V observations. Alternatively 
we consider the [N (N) / N (0)] uv value to be more reliable than the N(N) UV / N(R) V value 
given by equation (16), because for the UV line ratios the dependence on C(H/3) and T e is 
smaller and the slit size and position are the same. Similarly we consider [N (O) / N (H)] v to 
be more reliable than N(0)uv/N(H)v because the visual line ratios depend less on C(H[3) 
and T e than the UV to V line ratios, and the slit size and position are the same for the 
visual line intensities. 

To derive the same N/H value from equations (14) and (15) we need an ICF(N) twice 
as large as that provided by equation (14). A similar increase in the ICF(N) value was 
obtained for NGC 346, the brightest H II region in the SMC, by Relaho et al. (2002) based 
on CLOUDY ionization structure models. 

The O abundance was obtained from (Kingsburgh & Barlow 1994): 




( 



iV(He + ) + iV(He ++ ) 
iV(He + ) 



) 



2/3 



( 



iV(Q + ) + iV(0 ++ ) 
iV(H+) 



(17) 



where ICF{0) takes into account the fraction of +3 . 

The Ne abundance was obtained from (Peimbert & Costero 1969): 




(18) 



The S abundance was obtained from: 



JV(H) 



N(S) 



ICF(S) 



N(S + ) + iV(S ++ ) 
JV(H+) 



(19) 



where JCF(S) was estimated from the models by Garnett (1989) and amounts to 6.3 ± 1.5. 
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The CI abundance is given by: 

iVW - ICF(Cl) iV(H+) ' (20) 

where ICF(C\) was estimated from models computed with CLOUDY (Ferland 1996; Ferland 
et al. 1998) with similar O and He ionization structure and amounts to 1.05, a value very 
close to 1. the correction is due Cl +4 , the amount of Cl + for this object is negligible. 

Finally the Ar abundance is given by: 

^(Ar) Ar , iV(Ar ++ ) + iV(Ar+ 3 ) + N(Ar +4 ) 

= ICF{Al) n(^) ' (21) 

where ICF(Ar) = !/[!- 0+/0] = 1.02 (Liu et al. 2000), takes into account the Ar+ fraction. 



Discussion 



7.1. The A 4649 O II versus electron density diagram 



The relative intensity ratios of the O II recombination lines of multiplet 1 do not agree 
with the LTE computation predictions, this deviation was noted before for other objects by 
Esteban et al. (1999a); Esteban et al. (1999b); Peimbert (2003); and Tsamis et al. (2003). 
In particular Tsamis et al. (2003) discuss this problem extensively. 

To explain this deviation, we propose that the electron density is smaller than that 
needed to produce an LTE collisional redistribution of the fine structure energy levels. To 
test our proposal we produced Figure 2 that is based on the best O II observations in the 
literature, those presented in Table 8. We decided to plot 7(4649) over the sum of the 
intensities of all the lines of the multiplet /(sum), versus the density derived from optical 
forbidden line ratios, mainly the [CI III] lines that originate in the same regions where O is 
twice ionized. 

From Figure 2 we find an excellent correlation between /(4649)//(sum) and the for- 
bidden line electron density, iV e (FL). We propose to adjust the data with a curve given 
by 



7(4649) ~ 




7(4649)" 


/(sum) 


obs 


/(sum) 



+ 



"7(4649) 
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LTE 


" 7(4649)" 
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, iV e (crit) 
1 "T JV e (FL)_ 





(22) 



where [/(4649)//(sum)] cas is the value produced by recombination and cascade without col- 
lisional redistribution, [/(4649)//(sum)] L TE amounts to 0.397 Wiese et al. (1996) the value 
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predicted by collisional redistribution assuming that the populations of the upper energy 
levels of multiplet 1 are given by the statistical weights multiplied by the Einstein A val- 
ues, and iV e (crit) is the critical density defined by equation (22). From Figure 2 we find 
iV e (crit) = 1325±275 cm" 3 , and /(4649)//(sum)] cas = 0.105±0.015, therefore equation (22) 
can be written as: 



7(4649) 



292 

^ ai05+ i™/« (23) 



/(sum) 

It is beyond the scope of this paper to check this equation based on an atomic physics 
computation. 

Equation 22 has the following applications: a) when iV e (FL) is known it is possible to 
use this equation to obtain /(sum) from the observed /(4649) and from /(sum) to obtain 
the ++ abundance since /(sum) is independent the density, from /(4649) alone it is not 
possible to determine an accurate ++ abundance since this line is density dependent; b) 
from an observed /(4649)//(Sum) it is possible to determine N e in the ++ zone; c) if the 
forbidden line density of a nebula does not fall close to the value predicted with this equation 
it could mean that there are extreme inhomogeneities in density or chemical composition in 
this object, for example an object with high-density oxygen-rich knots embedded in a low 
density medium would fall to the left of the curve. 

From Table 8 it can be seen that the intensity of the lines that originate in the 3p 4 Dgy 2 
and 3p 4 D° / / 2 energy levels decreases with increasing local density, while the intensity of those 
lines that originate in the 3p 4 D°y 2 increases with increasing density. This is due to the effect 
of collisional redistribution that increases the population of the high statistical weight levels 
at the expense of the low statistical weight ones. Note that the intensity of the lines that 
originate from the 3p 4 D°^ 2 energy level are almost independent of the electron density. 

The results derived from Figure 2 imply that in these objects the [O III] AA 4363, 
4959, and 5007 lines and the O II recombination lines of multiplet 1 originate in regions of 
similar density, this is particularly the case for those objects with 7V e (FL) < 3000 cm -3 . In 
other words: we can not assume the presence of high density knots to explain the difference 
between the ++ abundances derived from the [O III] lines and the O II lines. Consequently 
we propose to drop the constant temperature assumption, and to explain the apparent 
difference in abundances as being due to the presence of temperature fluctuations. 

In general to reconcile the abundances derived from collisionally excited lines and re- 
combination lines it is necessary to assume the presence of large temperature variations. The 
presence of temperature variations in gaseous nebulae has been amply reviewed in the recent 
literature (e.g. Peimbert 2002; Esteban 2002; Liu 2002; Torres- Peimbert & Peimbert 2003). 
Even more recent results in favor of large temperature variations in the Orion nebula, the 30 
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Doradus nebula, and other planetary nebulae have been presented by Rubin et al. (2003); 
O'Dell et al. (2003); Peimbert (2003); Peimbert et al. (2003). 

Peimbert, Storey, & Torres-Peimbert (1993), by comparing the O II line intensities of 
multiplets 1, 2, and 10 for for the planetary nebula NGC 6572, have concluded that the O II 
lines are produced by recombination, and that the contribution of resonance fluorescence to 
the O II line intensities is not important. A similar conclusion was reached Grandi (1976) 
for the Orion nebula. 



7.2. The heliocentric velocities and the FWHM of the O II and [O III] lines 

The difference between the O(CEL) abundance and the O(RL) abundance indicates 
the presence of large temperature variations. There are two possible scenarios for these 
temperature variations: a) that the abundances are homogeneous along the line of sight and 
temperature fluctuations of ±23% are present, an explanation for these large temperature 
variations has to be sought (e.g. Torres-Peimbert & Peimbert 2003, and references therein), 
or b) that the abundances are not homogeneous along the line of sight and that there are 
high density knots with excess heavy element abundances embedded in a medium with 
lower heavy element abundances and higher electron temperature (e.g. Jacoby & Ford 1983; 
Torres-Peimbert, Peimbert & Pena 1990; Guerrero & Manchado 1996; Pequignot et al. 2002; 
Wesson et al. 2003). 

The second possibility is based on the idea that the central stars, of at least some 
planetary nebulae, firstly eject an envelope at relatively low velocity with almost normal 
heavy element abundances, and later on eject knots or clumps of material at higher velocities 
with an excess of heavy element abundances. In Abell 58 Guerrero & Manchado (1996) 
studied a hydrogen deficient high velocity knot blueshifted relative to the main body of the 
nebula by 120 km s -1 . while Wesson et al. (2003) studied three knots in Abell 30, where 
knot Jl is redshifted relative to knot J3 by about 45 km s _1 and knot J3 is redshifted relative 
to knot J4 by about 100 km s _1 . 

By assuming that chemical inhomogeneities are responsible for the temperature varia- 
tions in NGC 5307 we would expect most of the flux of the [O III] lines to originate in the 
lower density medium that expands at lower velocity, and a considerable fraction of the flux 
of the O II lines to originate from knots or clumps that move at higher velocities. Therefore 
if this were the case the radial velocity derived from the central wavelength of the [O III] 
lines could be different from that derived from the O II lines; furthermore the full width at 
half maximum, FWHM, of the [O III] lines would be smaller than that of the O II lines. 
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For NGC 5307 the heliocentric velocity of the AA 4363, 4959, 5007 [O III] lines is 30.1 
±0.5 km s -1 , while that of the O II lines of multiplet 1 is 33.1 ±4 km s -1 . The FWHM of 
the AA 4363, 4959, 5007 [O III] lines is 48.4 ±0.7 km s _1 , while that of the A 4649 O II line, 
the best observed line of multiplet 1, is 44.3 ±7 km s _1 . These results imply that there is 
no evidence in favor of high velocity O-rich clumps in our observations of NGC 5307. 



7.3. Comparison with other abundance determinations 

Table 7 presents a comparison between the abundance determinations of this paper and 
those by other authors. A discussion of the uncertainties in the determinations by other 
authors is beyond the scope of this paper, the errors of this work are presented in Tables 
6 and 9. The average JV(He)/iV(H) value derived by the other four groups amounts to - 
1.02 dex in excellent agreement with our value of -1.01 dex. The C abundance derived by 
us is based on the A 4267 C II recombination line, while those derived by other authors 
are based on UV lines, we consider our determination to be more accurate because it is 
almost independent of the adopted temperature while the abundance derived from the UV 
collisionally excited lines depends strongly on the adopted electron temperature (although 
the C II recombination abundances are still regarded as controversial by other authors). The 
N/H value has been already amply discussed and most of the difference among the various 
groups depends on the use of V or UV lines to determine the N abundance. Most of the 
difference among the O and Ne abundances is due to the adopted t 2 value. 

To discuss the effect of stellar evolution on the observed abundances in NGC 5307 we 
present in Table 9 the Orion nebula and solar abundances as representative of the present day 
interstellar medium and of the interstellar medium at the time the Sun was formed. Table 9 
presents the solar photospheric values for C, N, O, Ne, and Ar, and the solar abundances 
derived from meteoritic data for S, CI, and Fe. For the solar initial helium abundance the 
Y by Christensen-Dalsgaard (1998) was adopted, and not the photospheric one because, 
apparently, it has been affected by settling. Notice that the solar oxygen abundance used in 
this paper is the average value of the recent results by Allende-Prieto, Lambert, & Asplund 
(2001) and Holweger (2001) that indicate an O/H value of 8.71 dex, a value lower than that 
derived by Grevesse & Sauval (1998). We are also using the C/H value by Allende-Prieto, 
Lambert, & Asplund (2002), the N/H and Ne/H solar abundances were taken from Holweger 
(2001), while the S/H, Cl/H, and Ar/H come from Grevesse & Sauval (1998), the first two 
from the meteoritic data, while the last one from the photospheric data. 

The relative abundances of H, He, C, N, and O imply that NGC 5307 is a PN of Type 
II with a main sequence progenitor star in the 1.2 to 2.4 solar mass range (Peimbert 1978, 
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1990; Kingsburgh & Barlow 1994). 

In Table 9 we present the abundances of NGC 5307, the Orion nebula (representative of 
the interstellar medium today), and the Sun (expected to be similar to the interstellar abun- 
dances when the parental star of NGC 5307 was formed). By comparing these abundances 
we reach the following conclusions: a) the He abundance has not been significantly enriched 
by stellar evolution, b) the C abundance has decreased by about a factor of three, and c) the 
N abundance has increased by about a factor of three. These three results taken together 
indicate that the nebular material has only gone through part of the CNO cycle, converting 
C into N, but without producing significant amounts of helium. Evolutionary models of 
intermediate mass stars predict the N increase and the C decrease in their outer layers while 
the O and He abundances are almost unaffected (e. g. Torres- Peimbert & Peimbert 1971); 
N increases at the expense of C because the C/N ratio in the interstellar medium and in 
the Sun is considerably larger than the equilibrium value produced by the CNO cycle (e. g. 
Pagel 1997). 

Objects that show O-rich knots, like A30 and A58, show He/H excesses that are not 
present in NGC 5307 (Guerrero & Manchado 1996; Wesson et al. 2003). Consequently the 
chemical composition of NGC 5307 is also against the presence of O-rich knots in this object 

8. Conclusions 

NGC 5307 is a typical Type II PN with t 2 = 0.056 ±0.005. Its expansion velocity is low, 
about 15 km s -1 . Its chemical composition indicates that nuclear processing in the parental 
star of the ejected shell has not gone beyond the partial conversion of C into N. In particular 
it has a very similar He/H ratio to that present in the interstellar medium at the time its 
progenitor star formed 

Its chemical composition, together with the FWHM and the radial velocities of the O II 
lines and the [O III] lines, indicate that NGC 5307 does not have high density O-rich knots 
embedded in a low density O-normal medium. Consequently the difference between the O II 
and [O III] line intensities has to be due to a different cause than chemical inhomogeneities. 

The comparison of the FWHM of the O II and [O III] lines provides an excellent method 
to detect the presence O-rich knots, with a high electron density and a high velocity relative 
to the parental star, embedded in a low density O-normal medium. 

From the 1(4649) //(sum) versus iV e (Cl III) diagram we have found an equation that 
provides a good fit to most of the well observed gaseous nebulae. This fit to the data indicates 
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that most, if not all, of the O II recombination lines intensity originates in the same regions 
than the [O III] collisional lines intensity, and that high density O-rich knots are not the 
main cause responsible, in these objects, for the presence of temperature variations. 

We acknowledge the referee for a careful reading of the manuscript. MTR received par- 
tial support from FONDAP(15010003), a Guggenheim Fellowship and Fondecyt(1010404). 
MP received partial support from DGAPA UN AM (grant IN 114601). 
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Table 1. Journal of Observations 



Date A (A) Exp. Time (s) 



2002 March 12 3100-3880 3x90 

2002 March 12 3730-4990 3x180 

2002 March 12 4760-6840 3x90 

2002 March 12 6600-10360 3x180 



Table 2. Line Intensities 
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Table 2 — Continued 
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Table 2 — Continued 
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2 
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12 
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II 


1 
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25 
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4640.64 


N III 


2 


0.576 
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8 


4642.311 


4641.81 


II 


1 


0.187 
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15 
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N III 


2 








4647.945 


4647.42 


C III 


1 


0.049 


0.052 


25 


4649.767 


4649.13 


II 


1 
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0.215 


15 


4651.369 


4650.84 


II 


1 
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Table 2 — Continued 



A 


A 




Ion 


Id. 


pa, 


P 


err(%) 


4656 845 
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TfcW WW . 
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A 
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34 
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90 

^W 
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Np T? 
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W . 1 1 w 
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20 


4859 81 6 
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1 

1 w 
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33 
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11 1 


H4 
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1 00 

1 WW . w 
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1 

1 . w 
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Q3 


Hp T 
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48 
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w 
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Fl 
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10 
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[N II] 


F3 
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0.284 
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.67 


He I 
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2.0 
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.83 


[K IV] 


Fl 
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15 
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.70 


He II 
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0.079 


0.055 


25 


6301.241 
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30 


[0 I] 


Fl 


2.19 


1.49 


5 
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10 c 


[S III] 


F3 
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12 


6312.892 
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10 d 


[S III] 


F3 


1.20 
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7 


6364.728 


6363. 


.78 


[0 I] 


Fl 


0.865 


0.580 


8 


6394.571 


6393. 


6 


[Mn V] 




0.088 


0.059 


25 
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Table 2 — Continued 



A 


A 




Ion 


Id. 


pa, 


P 


/IT/ \ 

err(%) 
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Table 2 — Continued 
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He I 
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6/10 c 0.208 


0.098 


30 


9015.747 


9014.91 


H 
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Table 2 — Continued 



A 


Ao 


Ion 


Id. 




J b 


err(%) 


9069.106 


9068.9 C 


[S III] 


Fl 


6.78 


3.17 


6 


9070.260 


9068.9 d 


[S III] 


Fl 


12.1 


5.66 


4 


9211.235 


9210.28 


He I 


83 


0.224 


0.103 


30 


9229.907 


9229.02 


H I 


P9 


5.48 


2.52 


7 


9526.977 


9526.17 


He I 


6/8 c 


0.334 


0.150 


30 


9546.994 


9545.98 


H I 


P8 


6.02 


2.69 


7 


9560.951 


9560.3 


[Fe V] 




0.129 


0.058 


40 


10028.61 


10027.7 


He I 


6/7 c 


0.365 


0.157 


35 


10050.36 


10049.4 


H I 


P7 


9.27 


3.98 


7 



a Where F is the observed flux in units of 100.00 = 1.184 x 
10~ 12 erg cm -2 s _1 . 

b Where / is the reddening corrected flux observed flux, 
with C(H/3) = 0.59 dex, in units of 100.00 = 4.606 x 
10~ 12 erg cm -2 s _1 . 

c Blue component. 
d Red component. 

e See Pequignot & Baluteau (1988). 
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Table 3. Temperatures and Densities 



Lines T e (K) Lines iV e (cm 3 ) 



[S II] 


11230±1100 


[S II] 




[N II] 


11180±750 a 


[0 II] 


2120±240 


[0 II] 


15300±3000 b 


[CI III] 


2D4D+ 1350 


[S III] 


11020±360 


[Ar IV] 


2920±750 


[0 III] 


11800±100 






T(Bac) 


10700±2000 







a Corrected for recombination contribution to 
the A 5755 auroral line, see text. 

a Corrected for recombination contribution to 
the AA 7319,7320,7331, and 7332 auroral lines, see 
text. 
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Table 4. He + Ionic Abundance 



Line 


He+/H+ a 


3188 


7788±623 


3889 


9389±376 


4026 


8990±477 


4471 


8692±261 


4922 


8965±538 


5876 


9154±183 


6678 


9391±376 


7065 


9090±364 


Adopted 


9055±151 b 



a Given in units of 10 5 , for N e = 
2500, T 3 889 = 0.93, and t 2 = 0.031. 
Only the line intensity errors are 
presented. 

b In addition to the line intensity 
errors it also includes the effects of 
the uncertainties in N e = 2500 ± 
500, r 3889 = 0.93 ± 0.57, and t 2 = 
0.031 ±0.014. 
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Table 5. He, C, and O Ionic Abundances from Recombination Lines a 



Ion 


NGC 5307 


He+/H+ 


10.957±0.007 


He ++ /H+ 


9.86±0.01 


C++/H+ 


7.95±0.08 


0++/H+ 


8.77±0.03 


C+ 3 /H+ 


7.46±0.10 



a In units of 12+log(X m /H+). 
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Table 6. Ionic Abundances from Collisionally Excited Lines a 



Ion 


t 


2 = 0.00 


t 2 


= 0.056 


N+ 


6. 


24±0. 


.05 


6. 


,39±0 


,05 


0+ 


6. 


98±0. 


.08 


7. 


15±0. 


,08 


0++ 


8. 


48±0. 


.01 


8. 


77±0. 


.03 


Ne++ 


7. 


86±0. 


.03 


8. 


18±0. 


.04 


S+ 


5. 


03±0 


.05 


5. 


18±0. 


.05 


S++ 


6. 


20±0. 


.03 


6. 


42±0. 


.04 


C1++ 


4. 


51±0. 


.04 


4, 


79±0. 


.05 


C1+ 3 


4. 


,72±0. 


.03 


4. 


96±0. 


.04 


Ar++ 


5, 


62±0. 


.02 


5. 


86±0. 


.03 


Ar +3 


5. 


61±0. 


.03 


5, 


91±0. 


.04 


Ar+ 4 


3. 


53±0. 


.14 


3, 


78±0. 


.14 


K +3 


4. 


13±0. 


.12 


4. 


40±0. 


.12 



a In units of 12+log(X m /H+). 
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Table 7. NGC 5307 Abundance Determinations 51 



Element 


This paper 


This paper 


Tpp b 


FP C 


KB d 


M c 




t 2 = 0.000 


t 2 = 0.056 


t 2 = 0.035 


t 2 = 0.000 


t 2 = 0.000 


t 2 = 0.000 


He 


11.00 


10.99 


11.01 


10.94 


11.02 


10.96 


C 


8.09 


8.09 






8.21 


8.43 


N 


8.04 


8.33 


7.74: 


8.01: 


> 8.56 : 


8.16: 





8.51 


8.80 


8.64 


8.83 


8.63 


8.55 


Ne 


7.89 


8.21 


8.00 




7.89 


7.76 


S 


7.02 


7.24 




6.98 


6.47 


6.81 


CI 


5.00 


5.22 










Ar 


5.93 


6.20 






> 5.68 


6.14 



a In units of 12+ Log N(X)/N(K), gaseous content only. The C abundances derived in this 
paper are based on recombination lines, all the other heavy element abundances on collisionally 
excited lines. 

b Torres-Peimbert & Peimbert (1977). 
c de Freitas Pacheco et al. (1992). 
d Kingsburgh Barlow (1994). 
e Mal'kov (1998). 



Table 8. O II Line Intensity Ratios a 





NJ CI III) 


3p 4 D° , 

1 7/2 


3p 4 D° , 

1 5/2 


3p 4 D° /9 

1 3/2 


3p 4 D° 

^ 1/2 

( c ) 


T e (OIII) 




Object 


( cm -3 ) 


( b ) 


( c ) 


( d ) 


(K) 


Source 


NGC 604 


iooi^° 


0.134 


0.218 


0.392 


0.255 


8150 


Esteban et al. (2002) 


30 Doradus 


295±30 


0.158 


0.254 


0.372 


0.216 


10200 


Pcimbcrt (2003) 


30 Doradus 


480±100 


0.171 


0.288 


0.336 


0.204 


10100 


Tsamis et al. (2003) 


Ml 7 


1050±^ 


0.204 


0.319 


0.324 


0.153 


8200 


Tsamis et al. (2003) 


Ml-42 


169012°° 


0.327 


0.294 


0.246 


0.124 


9220 


Liu et al. (2001) 


NGC 5307 


20401^50 


0.291 


0.274 


0.272 


0.163 


11800 


This paper 


M8 




0.237 


0.263 


0.317 


0.183 


8050 


Esteban et al. (1999b) 


NGC 3576 


2700±1300 


0.240 


0.250 


0.241 


0.269 


8850 


Tsamis et al. (2003) 


NGC 7009 


3000±900 


0.333 


0.305 


0.218 


0.144 


10500 


Hyung & Aller (1995) 


NGC 6153 


3830±800 


0.345 


0.292 


0.259 


0.104 


9110 


Liu et al. (2000) 


M2-36 


4830±1000 


0.344 


0.287 


0.254 


0.115 


8380 


Liu et al. (2001) 


NGC 6543 


5000±1000 


0.332 


0.289 


0.230 


0.149 


8150 


Hyung et al. (2000) 


Orion 1 


5110±800 


0.337 


0.291 


0.251 


0.121 


8300 


Esteban et al. (1998) 


Orion 2 


6170±900 


0.329 


0.273 


0.244 


0.154 


8350 


Esteban et al. (1998) 


NGC 6572 


20000±4000 


0.359 


0.319 


0.232 


0.090 


10900 


Hyung et al. (1994a) 


NGC 5315 


500 00+ 50000 
ouuuu_ 25000 


0.390 


0.302 


0.214 


0.093 


8850 


Peimbert et al. (2003) 


IC 4997 


- 10 6 


0.352 


0.282 


0.216 


0.150 


13000 


Hyung et al. (1994b) 



limit f oo 0.397 0.300 0.202 0.101 ... Wiese et al. (1996) 



intensity ratios of lines belonging to multiplct 1, 3s 4 P - 3p 4 D°. The 8 lines of the multiplet are divided in 4 groups 
according to their upper energy level. The upper energy levels are shown in the column headers. 

b /(4649//(Sum). 

c /(4642 + 4676)//(Sum). 

d /(4639 + 4662 + 4696) //(Sum) . 

c /(4651 + 4674)//(Sum). 

'High density limit, in this case the line intensities are proportional to the Einstein A coefficients. 
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Table 9. NGC 5307, Orion, and Solar Total Abundances' 



Element 


NGC 5307 a 


Orion b 


Sun c 


12 + log He/H 


10.99±0.01 


10.99±0.01 


10.98±0.02 


12 + log O/H 


8.80±0.03 


8.72±0.06 


8.71±0.05 


log C/O 


-0.71±0.09 


-0.19±0.08 


-0.21±0.10 


log N/O 


-0.47±0.15 


-0.93±0.10 


-0.78±0.12 


log Ne/O 


-0.59±0.10 


-0.82±0.12 


-0.71±0.09 


log S/O 


-1.56±0.12 


-1.54±0.12 


-1.51±0.08 


log Cl/O 


-3.58±0.12 


-3.29±0.13 


-3.43±0.08 


log Ar/O 


-2.60±0.08 


-2.23±0.21 


-2.31±0.08 



a Gaseous abundances, values for t 2 = 0.056, obtained in 
this paper 

b Peimbert (1993); Esteban et al. (1998, 2002), values for 
t 2 = 0.024. The O and C abundances have been increased 
by 0.08 and 0.10 dex respectively to take into account the 
fractions of these elements trapped in dust grains. 

c Christensen-Dalsgaard (1998); Grevesse &; Sauval (1998); 
Allende-Prieto et al. (2001, 2002); Holweger (2001). 
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Fig. 1. — Section of the echelle spectrum including the O II multiplet 1 lines (observed 
fluxes) . 

Fig. 2. — Intensity ratio of the O II multiplet 1 line A 4649 relative to the sum of the 
intensities of the eight lines of that multiplet versus the local electron density derived from 
forbidden line ratios. The solid line represents the best fit to the data, see equation (23) and 
Table 8. 
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